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We demonstrate that a wide range of viable f{R) parameterizations (including the Hu & Saw- 
icki and the Starobinsky models) can be expressed as perturbations deviating from the ACDM 
Lagrangian. We constrain the deviation parameter b using a combination of geometrical and dy- 
namical observational probes. In particular, we perform a joint likelihood analysis of the recent 
Supernovae type la data, the Cosmic Microwave Background shift parameters, the Baryonic Acous- 
tic Oscillations and the growth rate data provided by the various galaxy surveys. This analysis 
provides constraints for the following parameters: the matter density flmo, the deviation from 
ACDM parameter b and the growth index 7(2). We parametrize the growth index 7(0) in three 
manners (constant, Taylor expansion around z = 0, and Taylor expansion around the scale factor). 
We point out the numerical difficulty for solving the generalized f{R) Friedman equation at high 
redshifts due to stiffness of the resulting ordinary differential equation. We resolve this problem by 
constructing an efficient analytical perturbative method in the deviation parameter b. We demon- 
strate that this method is highly accurate, by comparing the resulting analytical expressions for 
the Hubble parameter, with the numerical solutions at low and intermediate redshifts. Surprisingly, 
despite of its perturbative nature, the accuracy of the method persists even for values of b that are 
of 0(1). 

PACS numbers: 98.80.-k, 98.80.Bp, 98.65.Dx, 95.35.-|-d, 95.36. +x 



1. INTRODUCTION 

A variety of cosmological studies have converged to 
a cosmic expansion history involving a spatially flat ge- 
ometry and a cosmic dark sector formed by cold dark 
matter and some sort of dark energy, endowed with large 
negative pressure, in order to explain the observed ac- 
celerating expansion of the Universe [ll-Q. In this frame- 
work, the absence of a fundamental physical theory, re- 
garding the mechanism inducing the cosmic acceleration, 
has given rise to a plethora of alternative cosmological 
scenarios. Modified gravity models act as an important 
alternative to the scalar-field dark energy models, since 
they provide an efficient way towards explaining the ac- 
celerated expansion of the universe, under a modification 
of the nature of gravity. Such an approach is an attempt 
to evade the coincidence and cosmological constant prob- 
lems of the standard ACDM model. 

Particular attention over the last decades has been paid 
on f{R) gravity theories [l3|- In this scenario of non- 
standard gravity, one modifies the Einstein-Hilbert ac- 
tion with a general function f{R) of the Ricci scalar R. 



The /(i?) approach is a relative simple but fundamen- 
tal tool used to explain the accelerated expansion of the 
universe. A pioneering approach was proposed long ago, 
where f{R) = R + mR^ Later on, the f{R) mod- 

els were further explored from different points of view in 
[l^ - [l3 | and a large number of functional forms of f{R) 
gravity is currently available in the literature. It is in- 
teresting to mention that subsequent investigations [l3| 
confirmed that 1/R gravity is an unacceptable model be- 
cause it fails to reproduce the correct cosmic expansion 
in the matter era. Of course, there are many other pos- 
sibilities to explain the present accelerating stage. In- 
deed, in the literature one can find a variety of modified 
gravity models (for reviews see Refs. flE'l) which include 
the braneworld Dvali, Gabadadze and Porrati (hereafter 
DGP; [16]) model, Finsler-Randers gravity [I3, scalar- 
tensor theories [1^ and Gauss-Bonnet gravity jl9| . 

The construction of observationally viable f{R) the- 
ories has proved to be challenging because it has been 
shown [14| that most of these models do not predict a 
matter era in the cosmic expansion history. Neverthe- 
less observationally viable f{R) models have been con- 
structed and two such examples are the following: 

1. The Hu & Sawicki model [2lj with 
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where Ci, C2 are free parameters, ~ ^moH^ is of 
the order of the Ricci scalar Rq at the present time, 
Hq is the Hubble constant, flmo is the dimension- 
less matter density parameter at the present time, 
m and n are positive constants. 



2. The Starobinsky model [2g with 



f{R) = R-cim^ 1 - (1 + i?V 



(1.2) 



These models were originally advertised as models that 
do not contain the cosmological constant as part of f{R) 
being distinct from the ACDM form f{R) = R ~ 2A 
(where A is the cosmological constant). However, it is 
straightforward to show that both the Hu & Sawicki 
and the Starobinsky models may be written in terms of 
ACDM modified by a distortion function y{R, b) that de- 
pends on a deviation parameter b as: 



fiR)^R-2Ay{R,b) 



where |2C 



yiR,b) = i 



1 



1 + {R/ib A)" 



(1.3) 



(1.4) 



for the Hu & Sawicki model with A = ^^^-^ and b = 

2C2 

2cl 



-, while 



2/(i?,6) = l- 



(l + (^)T 



(1.5) 



for the Starobinsky model where A = and b — 

Notice that in both cases the following two limits exist 
for n > 0: 



lim/(i?) = R-2A 
lim f{R) = R 



(1.6) 



and therefore both models reduce to ACDM for b ^ 0. 
Notice that both the Hu & Sawicki and the Starobinsky 
models effectively include the cosmological constant even 
though they were advertised as being free from a cosmo- 
logical constant in the original papers [2l|, . In fact 
by proper choices of the function y{R, b) it is possible to 
construct infinite viable f{R) models which however will 
always include ACDM limiting case for 6—^0. 

Thus, an important question that arises is the follow- 
ing: 'What is the range of the deviation parameter b that 
is consistent with cosmological observations?' This is 
the main question addressed in the present study. Since 
ACDM is consistent with observations we anticipate that 
the value 6 = is within the acceptable range of b val- 
ues. Thus, the interesting part of the question is 'What is 
the maximum allowed value of b at eg the 2a confidence 
level?' 



In order to address this question we solve the back- 
ground modified Friedman equation assuming flatness 
and obtain the Hubble parameter H{nmo, b] z). This in- 
volves the numerical solution of a stiff ordinary differ- 
ential equation (ODE), second order in iJ, with initial 
conditions at high z that correspond to ACDM . The 
numerical solution of this stiff ODE at redshifts higher 
than z ~ 300 is quite challenging. However, we have 
developed an efficient analytical perturbative expansion 
in b to solve it. This expansion leads to an analytic ex- 
pression for H{njno, b; z) to all orders in b. We thus use 
geometric probes (type la surenovae, CMB shift param- 
eter and Baryon Acoustic Oscillation data) to constrain 
the parameters flrnQ,b that appear in the expression of 
H{nrno,b;z). 

In addition to geometric observations that probe di- 
rectly the cosmic metric, dynamical probes play a crucial 
role in constraining cosmological models. The growth in- 
dex, 7, could provide an efficient way to discriminate 
between modified gravity models and scalar field dark 
energy (hereafter DE) models which adhere to general 
relativity. The accurate determination of the growth 
index is considered one of the most fundamental tasks 
for Observational Cosmology. Its importance steams 
from the fact that there is only a weak dependence of 
7 on the equation of state parameter u;(z), as has been 
found by Linder & Cahn [23j , which implies that one can 
separate the background expansion history, H{z), con- 
strained by geometric probes (SNIa, BAO, CMB), from 
the fluctuation growth history, given by 7. For a con- 
stant DE equation of state w, it was theoretically shown 
that for DE models within general relativity the growth 
index 7 is well approximated by 7 ~ ^g^Zs^ (^^e Refs. 
nil,!!!,!!!],!!!!), which boHs down to Z 6/11 for the 
ACDM cosmology w{z) = —1. Notice, that in the case of 
the braneworld model of Dvali, Gabadadze & Porrati [l^ 
we have 7 « 11/16 (see also Refs. pi l27l-[29l|). while for 



0.415 - 0.21z for 
3l|). Recently, Basi- 
9/14 for the Finsler- 



some f(R) gravity models we have 
various parameter values (see [13, 
lakos & Stavrinos [11] found 7 « 
Randers cosmology. 

Observationally, indirect methods to measure 7 have 
also been developed (mostly using a constant 7) , based ei- 
ther on the observed growth rate of clustering ^26i,i27i. 
|40| providing a wide range of 7 values 7 ~ (0.58 — 0.67), 
or on massive galaxy clusters (Vikhlinin et al. |4l| 
and Rapetti et al. [13] )■ The latter study provides 
7 = 0.42tQ;^g. An alternative method for measuring 7 
involves weak gravitational lensing Gaztanaga et 

al. [l^l performed a cross-correlation analysis between 
probes of weak gravitational lensing and redshift space 
distortions and found no evidence for deviations from 
general relativity. Also, Basilakos & Pouri [13| and Hud- 
son & TurnbuU [3811 imposed constraints on the growth in- 
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dex using the combination parameter F{z)as{z) ^, of the 
growth rate of structure F{z) muhiphed by the redshift- 
dependent rms fluctuations of the hnear density field, 
as{z). The above authors found 7 = 0.602 ± 0.05 fST^ 
7 = 0.619 ±0.05 [Hi while Basilakos & Pouri 37] showed 
that the current growth data can not accommodate the 
Dvali, Gabadadze & Porrati ^16] gravity model. 

In order to impose constraints on the viable f{R) 
models discussed above, we use in addition to geomet- 
ric probes, the recent growth rate data as collected by 
Nesseris & Garcia-Belhdo [i^, Hudson & Turnbull fssj 
and Beutler et al. [46j . 

The plan of the paper is as follows. Initially in sec- 
tion II, we briefly discuss the background cosmological 
equations. The basic features of the growth index are 
presented in section III, where we extend the original Po- 
larski & Gannouji method [l^l for a general family of 7(2) 
parametrizations as well as f{R) cosmological models. In 
section IV, a joint statistical analysis based on the Union 
2.1 set of type la supernovae (SNIa the observed 

Baryonic Acoustic Oscillations (BAOs \^\), the shift pa- 
rameter of the Cosmic Microwave Background (CMB [9|), 
and the observed linear growth rate of clustering, mea- 
sured mainly from the PSCz, 2dF, VVDS, SDSS, 6dF, 
2MASS, BOSS and WiggleZ redshifts catalogs, is used 
to constraint the growth index model free parameters. 
Finally, we summarize our main conclusions in section 
V. 



2. THE BACKGROUND EVOLUTION 

First of all we start with the assumption that the uni- 
verse is a self-gravita ting fluid described by a modified 
gravity namely f{R) and endowed with a spatially 
flat homogeneous and isotropic geometry. In addition, we 
also consider that the universe is filled by non-relativistic 
matter and radiation. The modified Einstein-Hilbert ac- 
tion reads: 



where R is the Ricci scalar, fn = df /dR, is the Ein- 
stein tensor and Tjf is the energy-momentum tensor of 
matter. Modeling the expanding universe as a perfect 
fluid that includes radiation and cold dark matter with 
4- velocity C/^, we have Tif = -Pg^; + {p + P)U^'U^, 
where p — /?,„ -f pr and P = Pm + Pr are the total en- 
ergy density and pressure of the cosmic fluid respectively. 
Note that pm is the matter density, pr denotes the den- 
sity of the radiation and = 0, = Pr /3 are the corre- 
sponding pressures. Assuming negligible interaction be- 
tween non-relativistic matter and radiation the Bianchi 
identity v'^ T^i, = (which insures the covariance of the 
theory) leads to the matter/radiation conservation laws: 



Pn 



iHpjn = Pr 



iHpr- = 

and Pr 



(2.3) 



the solutions of which are Pm = PmoQ and. p^ = Proo- 
Note that the over-dot denotes derivative with respect to 
the cosmic time t, a{t) is the scale factor and H = a/a is 
the Hubble parameter. 

Now, in the context of a flat FLRW metric with Carte- 
sian coordinates 

ds'^ = -df + a^{t){dx^ + dy'^ + dz^) (2.4) 

the Einstein's tensor components are given by: 



Si' 2H + 3H' 



(2.5) 



Inserting Eqs. (|2.5p into the modified Einstein's field 
equations (|2.2p . for comoving observers, we derive the 
modified Friedmann's equations 

3/i^i^' - + iHfuRR = {pm + Pr) (2.6) 



- = e [prn + (4/3)p,] + fu. - HfR (2.7) 

where R = aHdR/da and fffR = d'^f/dR^. Also, the 
contraction of the Ricci tensor provides the Ricci scalar 



S 



d'^x. 



1 
2k 



J{R}+Cra+Cr 



(2.1) 



where Cm is the Lagrangian of matter, Cr is the La- 
grangian of radiation and k'^ — 8ttG. Now varying the 
action with respect to the metric^ we arrive at 



2nfR -if- Rfu) 



(2.2) 



^ Here the capital F{a) denotes the growth rate of structure. We 
follow the latter notation in order to avoid confusion with the 

fiR)- 

^ We use the metric i.e. the Hilbert variational approach. 



R = gf^'R^, = 6 



6{2H^ + H). (2.8) 



Of course, if we consider f{R) = R then the field equa- 
tions ()2.2|) boil down to the nominal Einstein's equations 
a solution of which is the Einstein de Sitter model. On 
the other hand, the concordance A cosmology is fully re- 
covered for f{R) = R — 2A. We would like to stress 
here that within the context of the metric formalism the 
above f{R) cosmological models must obey simultane- 
ously some strong conditions (for an overall discussion see 
15]). Briefly these are: (i) /jj > for i? > i?o > 0, where 
Rq is the Ricci scalar at the present time. If the final at- 
tractor is a de Sitter point we need to have /r > for 
R> Ri > 0, where i?i is the Ricci scalar at the de Sitter 
point, (ii) Jrr > for i? > i?o > 0, (iii) /(i?) « i? - 2A 
for i? > i?o and finally (iv) < ^jf^{r) < 1 at 

r - -Mil - _9 
r- f - ^- 
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Finally, from the current analysis it becomes clear that 
unlike the standard Friedmann equations in Einstein's 
general relativity the modified equations of motion ()2.6|) 
and (12.71) are complicated and thus it is difficult to solve 
them analytically. Below, we are going to compare the 
f{R) results with those of the concordance ACDM model. 
This can help us to understand better the theoretical 
basis of the current f{R) models as well as the variants 
from general relativity. For practical reasons (see below) , 
we derive the the effective ("geometrical") dark energy 
EoS parameter in terms of E{a) — H{a)/HQ (see [Tsllsgl 
and references therein) 



where 



n,„(a) 



-1 



2 ^ rflnE 

3 da 



1 - Qr,i{a) 



E^{a) 



(2.9) 



(2.10) 



Differentiating the latter and utilizing Eq. (|2.9|) we find 
that 



da 



-w{a)nm{a) [1 - fim(a)] 



(2.11) 



In the case of the traditional ACDM cosmology f{R) — 
R — 2A, the corresponding dark energy EoS parameter is 
strictly equal to —1 and the normalized Hubble function 
in the matter era is given by 



1/2 



(2.12) 



A. The f{R) functional forms 

In order to solve numerically the modified Friedmann 
equation (|2.6|) we need to know apriori the functional 
form of f{R). Due to the absence of a physically well- 
motivated functional form for the f{R) parameter, there 
are many theoretical propositions in the literature. In 
this article for the background we use different reference 
expansion models namely flat ACDM and f{R) respec- 
tively. Bellow we briefly present the two most popular 
f{R) models whose free parameters can be constrained 
from the current cosmological data. 

Firstly, we use the Hu & Sawicki {21'] model (here- 
after /iCDM) as expressed by equation (|l.ip . Using the 
constrains provided by the violations of weak and strong 
equivalence principle, Capozziello & Tsujikawa [51] found 
that n > 0.9. On the other hand it has been proposed 



by [54I that n is an integer number, so for simplicity in 



our work we have set n = 1. In Ref. [21| the parameters 
(ci,C2) were related to ftmo, ^ro and the first derivative 
of f{R) at the present epoch fnQ in order to ensure the 
expansion history is close to that of ACDM. Specifically, 



we have 
£1 

C2 

Irq 



6 



(1 — D,rO — i^mo) 



mO 



nci 



12 i2ri^Q 



mO 



n 



mO 



The first two derivatives of Eq. (|l.ll) with respect to R 



R 


^2 ( ^^-^2 J 


n 

4 


- 1 


2 




R 


C2 


(AT 


+ 1 


2 



Irr 





C2 


[n + l) 


(4 

V 


-) -n + l 


i?2 


C2 




f 1 


3 



(2.13) 



(2.14) 



As discussed in the Introduction, the Lagrangian of 
Eq. (jl.ip can also be written as 



/(i?) = R 



m?ci 



m?ci/c2 



= i? - 2A 1 - 



= R- 



C2 l+C2(i?/m2)" 

1 



1 + {R/{b A)" 



2A 



where A = ^^^^^ and b 



1-1/" 



(2.15) 



In this form it is clear 



that the HS model can be arbitrarily close to ACDM, 
depending on the parameters b and n. 

We now consider the Starobinsky [l^] model (hereafter 
/2CDM) as expressed by equation (|1.2p : As in the Hu & 
Sawicki [2l| model we choose m? ~ flmoHo, while [5l| 
also showed that n > 0.9. In this case the fu and fnu 
derivatives are given by 



2nR 



fR^l 



Cl 



(2.16) 



Irr — — 



1 + ^j [to-* _ (1 + 2n)i?2] Cl 



(m4 + ii2y 

(2.17) 

In order to ensure that the expansion history of this 
model is close to that of ACDM we need to match the ci 
constant to A, ie — ci — — 2A = —6(1 — flmo — ^rojHg 
or 



Cl = 



6(1 — ilmO — ^ro) 



(2.18) 



As discussed in the Introduction, the Lagrangian of 
Eq. (|1.2p can also be written as 



f{R) = i? - 2A 



1 



(2.19) 
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where A = and b = —. In this form it is clear that 
that this model can also be arbitrarily close to ACDM, 
depending on the parameters b and n. Thus, the param- 
eter b determines how close the model is to ACDM. 

It is interesting to mention that the above f{R) models 
satisfy all the strong conditions (see section II) and thus 
they provide predictions which are similar to those of 
the usual (scalar field) DE models, as far as the cosmic 
history (presence of the matter era, stability of cosmo- 
logical perturbations, stability of the late de Sitter point 
etc.) is concerned. Also, we will restrict our present nu- 
merical solutions to the choice Hq = 70.4Km/s/Mpc and 
(Tg = 0.8. For example, in this case the modified Fried- 
mann equations for the f{R) models contain two free 
parameters, namely (r2„io,6) which can be constrained 
from the current cosmological data. 



B. Analytic approximations 

In this subsection we present a novel approximation 
scheme for the solution of the modified Friedmann equa- 
tion (j2.6p and we explicitly apply it to the two widely 
used models Eqs. ((2?T5l) and (|239)) . 

In particular we may write eq. (j2.6p as 

-IrH^N) + (l^^oe-'^ + ^roe-^^) + ^{JbR ~ /) 
= fRBH^{N)R'{N) , (2.20) 

where the prime denotes differentiation with respect to 
N and R{N) is given by eq. (|2.8p . Using now eq. (II. 3p 
for specific f{R) models, the above ODE (and its solution 
H{N)) may be expanded around ACDM with respect to 
the deviation parameter b. 

Since we are interested in testing deviations from the 
ACDM model, we find it useful to perform a series ex- 
pansion of the solution of the ODE (|2.20p around b — 
as 

M 

H'^iN) = Hl{N) + b'SHf{N) (2.21) 
1=1 

where 

— T72~^ ~ f^mOe '^'^ +ilrOe + {I ~ flmO~ ^ro) (2.22) 

and M is the number of terms we keep before truncat- 
ing the series. Usually keeping only the two first non-zero 
terms is more that enough to have excellent agreement of 
better than 0.001% at all redshifts with the numerical so- 
lution for realistic values of the parameter b £ [0.001, 0.5]. 

By expanding Eq. (I2.20p with (|2.2ip to any given or- 
der in b we can find analytical solutions for the Hubble 
expansion rate. It is easy to show that for the HS model 
and for n — I the first two terms of the expansion are 
the following: 



HjjsW ^ Hl{N) + b 5Hl{N)W 5Hl{N) + ... (2.23) 

where 5Hl{N) and dH^iN) are given by Eqs. (lA.lP and 
(|A.2p respectively. For the Starobinsky model for n = 1 
we have 



Hl,,,{N) = Hl{N) + b' SHliN) + h^ 5Hl{N) + ... 

(2.24) 

where 6H2{N) and 6H^{N) are given by Eqs. (TOI) and 
(|A.4p respectively. 

Obviously, similar expressions can be obtained for any 
f{R) model and up to any order provided that for & — >■ 
we obtain ACDM . We stress that the expressions for 
6Hf{N) are algebraic up to all orders, something that 
makes this method very useful and fast compared to solv- 
ing the differential equation numerically. Furthermore, 
this method avoids another problem of the numerical in- 
tegration, namely that at very high redshifts the ODE of 
Eq. (|2.6p is quite stiff, thus making the integration im- 
possible with standard methods. This makes the numeric 
solution quite time consuming and possibly unreliable, 
something which as we will show is not a problem for our 
analytic approximation. 

In what follows we will test the validity of this approx- 
imation in both cases. In order to do this we compared 
the predictions of the analytical solutions of Eqs. (I2.23P 
and (|2.24p to the numerical solution in each case for a 
large variety of values for the parameter b. In particular, 
we estimated the average percent deviation between the 
approximations and the numerical solution, defined as: 

(error(6)) = (lOQ.fl-pEEI^d^)) (2.25) 

where the average is taken over redshifts in the range 
z G [0, 30]. The reason for averaging is that most of the 
data we will use involve distance scales, like the luminos- 
ity distance, that are integrals of the Hubble parameter. 
We have also kept z below 200 since the numerical ODE 
solver of Mathematica is unable to go to larger redshifts 
due to the stiffness of the ODE. 

We show the results of the average error for a large 
variety of values of the parameter 6 e [0.01, 2] in Fig. [1] 
Clearly, the approximation behaves exceptionally espe- 
cially within the ranges of interest, ie the vertical dashed 
lines. These regions correspond, as we will see in a later 
section, to the best fit values of b. The horizontal dotted 
line indicates an error of 10~^%. Finally, we see that our 
approximation is on average in agreement to better than 
^ 0.01% for realistic parameters, ie 6 ~ 0(1), of the f{R) 
models. In a forthcoming paper we attempt to investi- 
gate the validity of our method against all the available 
f{R) gravity models. 
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FIG. 1: The results of the average error in the redshift range z G [0,30] for a large variety of values of the parameter 
h G [0.01,2]. Clearly, the approximation behaves exceptionally especially within the ranges of interest, ie the vertical dashed 
lines. The horizontal dotted line indicates an error of 10~^%. 



3. THE EVOLUTION OF THE LINEAR 
GROWTH FACTOR 

In this section we concentrate on the sub-horizon scales 
in which the DE component is expected to be homoge- 
neous and thus we can use perturbations only on the 
matter component of the cosmic fluid [53]. Therefore, 
the evolution equation of the matter fluctuations, for 
cosmological models where the DE fluid has a vanishing 
anisotropic stress and the matter fluid is not coupled to 
other species (see [13, [il,!!!,^!,!!^,!!,!!), is given 
by: 



(3.1) 



where is the matter density and G'off(i) = GjqQ(t)^ 
with Gat denoting Newton's gravitational constant. 

For those cosmological models which adhere to general 
relativity, = 1, Geff = Gat], the above equation 

reduces to the usual time evolution equation for the mass 
density contrast [60| , while in the case of modified gravity 
models (see [il,^, l^,!!!), we have Geff|^ Gn (or 
Q ^ 1). Indeed it has been shown (see Refs. ^3,113) that 
in the case of /(i?) models the quantity Q is a function 
of the scale factor and of the wave-number k 



(5(a, k) 



1 



1 + 4^ 



1 + 3^ 



■ fn 



(3.2) 



1/A 



We restrict our analysis to the choice of k 
0.1/iMpc-i or A = lO/i-^Mpc (see also lU). 

In this context, 5m{t) oc D{t), where D{t) is the lin- 
ear growing mode (usually scaled to unity at the present 
time). Of course, solving Eg. p.ip for the concordance A 
cosmology, we derive the well known perturbation growth 
factor (see [60l|l: 



Da{z) 



4-c>o 



(1 + u)du 
Eliu) 



(3.3) 



In this work we use the above equation normalized to 
unity at the present time. 

Since in most of the cases Eq. (13. ip does not yield an- 
alytical solutions, it is common in this kind of studies to 
provide an efficient parametrization of the matter pertur- 
bations which is based on the growth rate of clustering 

(eol 



F{a) 



d In 5„ 
din a 



(3.4) 



where 7 is the growth index (see Refs. |23l-l26ll54j) which 
plays a key role in cosmological studies as we described 
in the introduction. 



A. The generalized growth index parametrization 



Inserting the first equality of Eq. 
and using simultaneously Eq 
derive after some algebra, that 



into Eq. ((XT|) 



23 and ^^i?^, 



we 



dF 3 
a— + F'^+X{a)F=-VLra(a)Q{a), (3.5) 
da 2 



with 



X{a)^]^-^w{a) [l-a„(a)] 



(3.6) 



we 



where in order to evaluate the final form of Eq. 
have used Eq. (PTT|) . 

Now, we consider that the growth index varies with 
cosmic time. Transforming equation p.5p from a to red- 
shift [-^ = -(1 + z)-^£] and utilizing Eqs.dSH), (PTT|) 
we simply derive the evolution equation of the growth 
index 7 — 7(2) (see also fi^). Indeed this is given by: 



-(l + z)7'ln(r!. 



r!7„ + 3u,(l-f^„0(7-^) + ^ 

(3.7) 



7 



where prime denotes derivative with respect to redshift. 
At the present time the above equation becomes 



-7'(0)ln(17„o) + ^Zio^ + 3wo{l - f^™o)[7(0) 



2^ + 2 
(3.8) 



where Qq — Q{z = 0, k) and Wq — w{z — 0). 

In this work we phenomenologically parametrize 7(2:) 
by the fohowing general relation (see I37|) 



7(^) = lo+liviz) ■ 



(3.9) 



Obviously, the above equation can be viewed as a first 
order Taylor expansion around some cosmological quan- 
tity such as a{z), z and ^^{z). We would like to stress 
that for those y{z) functions which satisfy j/(0) = [or 
7(0) = 7o] one can write the parameter 71 in terms of 79. 
In this case [7'(0) — 7i2/'(0)], using Eg. p.Sp we obtain 



71 



0^»o + ^Mlo ~ i)(l - a„o) - iQof^^o"" + 



y'(0)lnf2„o 

(3.10) 

Let us now briefly present various forms of 7(2;), Vz. 

• Constant growth index (hereafter Tq model): Here 
we set 71 strictly equal to zero, thus 7 = 70. 

• Expansion around z = (see [l^]; hereafter Fi 
model): In this case we have y{z) = z. Note how- 
ever, that this parametrization is valid at relatively 
low redshifts < z < 0.5. In the statistical anal- 
ysis presented below we utilize a constant growth 
index, namely 7 = 70 + O.571 for z > 0.5. 



• Expansion around a — 1 (|62h64|: hereafter r2 
model): Here the function y becomes y{z) = 1 — 
^(■2) = Obviously, at large redshifts z 1 we 
get 7oo 70 + 71- 

To conclude, for the Fi and r2 parametrizations one can 
easily show that y{0) — and y'{0) = 1, respectively. 
As an example, for the case of the ACDM cosmology 
with 70 ~ 6/11 and flmo = 0.273, Eg. ljXTU)) provides 
71 ~ —0.0478. In addition, based on the Starobinsky's 
f{R) model with (f7™o,7o) (0.273,0.415) Gannouji et 
al. tSfl] found 71 = -0.21. 



4. OBSERVATIONAL CONSTRAINTS 

In the following we briefly present some details of the 
statistical method and on the observational sample that 
we adopt in order to constrain the free parameters of the 
growth index, presented in the previous section. 

A. The Growth data 



and WiggleZ galaxy surveys, for which their combina- 
tion parameter of the growth rate of structure, F[z)^ 
and the redshift-dependent rms fluctuations of the linear 
density field, (J%{z)^ is available as a function of redshift, 
F{z)a^{z). The Fa^ = /erg estimator is almost a model- 
independent way of expressing the observed growth his- 
tory of the universe (see |66[). Indeed the observed 
growth rate of structure {Fobs = f3B) is derived from the 
redshift space distortion parameter /3(z) and the linear 
bias B. Observationally, using the anisotropy of the cor- 
relation function one can estimate the /3(z) parameter. 
On the other hand, the linear bias factor can be defined 
as the ratio of the variances of the tracer (galaxies, QSOs 
etc) and underlying mass density fields, smoothed at 
8h~^ Mpc B{z) = as,tr{z)/o'8{z), where ag,,triz) is mea- 
sured directly from the sample. Combining the above def- 
initions we arrive at fag = Fag = j3ag,_tr- We would like 
to point that the different cosmologies (including those of 
modified gravity) enter only weakly in the observational 
determination of j3{z) (and thus of /erg), through the 
definition of distances. In Table I we quote the precise 
numerical values of the data points with the correspond- 
ing errors and references. 



B. The overall Likelihood analysis 

In order to constrain the cosmological parameters and 
the growth index of the /(i?) models one needs to per- 
form a joint likelihood analysis, involving the cosmic ex- 
pansion data such as SNIa and BAO, CMB shift param- 
eter together with the growth data. Up to now, due to 
the large errors of the growth data with respect to the 
cosmic expansion data, various authors preferred to con- 
strain first r^mo using SNIa/BAO/CMB and then to use 
the growth data alone. Of course, armed with the recent 
high quality growth data it would be worthwhile to si- 
multaneously constrain (rimO: 7). In particular, we use 
the Union 2.1 set of 580 SNIa of Suzuki et al. (H] ^, and 
the observed BAOs. The BAO data are given in terms of 
the parameter d^{z) = ^-^^^^f^, where lBAo{zdrag) is 
the BAO scale at the drag redshift, assumed known from 
CMB measurements, and Sal 



Dv(z) = 



{\ + zfDA{zf 



H{z) 



1/3 



(4.1) 



is the usual volume distance. 

In this analysis we use the 6dF, the SDSS and WiggleZ 
BAO data shown in Table [III The WiggleZ collaboration 
[7^ has measured the baryon acoustic scale at three dif- 
ferent redshifts, complementing previous data at lower 
redshift obtained by SDSS and 6DFGS 



The growth data that we utilize in this article based 
on the PSCz, 2dF, VVDS, SDSS, 6dF, 2MASS, BOSS 



^ The Snia data can be found in http;/ /supernova.lbl. gov/Union/ 
and in El 
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TABLE I: Summary of the observed growth rate and references. 



Index 


z 


growth rate (/cr«„jj 


Refs. 


1 


0.02 


0.360 ± 0.040 


[38] 


2 


0.067 


0.423 ± 0.055 


[46] 


3 


0.17 


0.510 ±0.060 


[65, 66] 


4 


0.35 


0.440 ± 0.050 


[66, 67] 


5 


0.77 


0.490 ±0.180 


[34, 66] 


6 


0.25 


0.351 ±0.058 


[68] 


7 


0.37 


0.460 ± 0.038 


[68] 


8 


0.22 


0.420 ± 0.070 


[69] 


9 


0.41 


0.450 ± 0.040 


[69] 


10 


0.60 


0.430 ± 0.040 


[69] 


11 


0.78 


0.380 ± 0.040 


[69] 


12 


0.57 


0.427 ± 0.066 


[70] 


13 


0.30 


0.407 ± 0.055 


m 


14 


0.40 


0.419 ± 0.041 


m 


15 


0.50 


0.427 ± 0.043 


m 


16 


0.60 


0.433 ± 0.067 


m 



TABLE IL The BAO data used in this analysis. The first six 
da ta p oints are volume averaged and correspond to Table 3 
of [T^l. Their inverse covariance Matrix is given by (|4.3p . 





6dF 


SDSS 


WiggleZ 


z 


0.106 


0.2 0.35 


0.44 0.6 0.73 


dz 


0.336 


0.1905 0.1097 


0.0916 0.0726 0.0592 


Adz 


0.015 


0.0061 0.0036 


0.0071 0.0034 0.0032 



correspond to the uncertainty of the Hubble parameter 
Hq today, i.e. Hq = 100 /ikmsec^^ Mpc"^). 

For a flat prior, the 9-year WMAP data (WMAP9) 
measured best-flt values are 1741 



CMB 




302.40 
1.7246 
1090.88 



(4.5) 



The corresponding inverse covariance matrix is [74 



The chi-square is given by 

xIao = TM - d{z,)]C-/[d, - d{z,)], (4.2) 



'CMB 



3.182 
18.253 
-1.429 



18.253 
11887.879 
-193.808 



-1.429 
-193.808 
4.556 



C. 



where the indices i,j are in growing order in z, as in 
Table ini For the first six points, C,^^ was obtained from 
the covariance data in . 7^ in terms of dz : 

0. 
0. 
0. 

10586 
-59907 
125536 
(4.3) 

The positions of CMB acoustic peaks are affected by 
the expansion history of the Universe from the decoupling 
epoch to today. In order to quantify the shift of acous- 
tic peaks we use the data points (^q, 7?., Zcmb) of Ref. [tS ] 
(WMAP9), where la and 7^ are two CMB shift parame- 
ters and Zcmb is the redshift at decoupling. 

For the redshift Zcmb there is a fitting formula by Hu 
and Sugiyama 75]: 



We thus define 

-X"cMB 



la - 302.40 
7e- 1.7246 
Zcmb - 1090.8 



/ 4444 


0. 


0. 


0. 


0. 


0. 


30318 


-17312 


0. 


0. 


0. 


-17312 


87046 


0. 


0. 


0. 


0. 


0. 


23857 


-22747 


0. 


0. 


0. 


-22747 


128729 


V 0. 


0. 


0. 


10586 


-59907 



Vmd construct the contribution of CMB to as 



XCMB 



"CMB ^-'CMB ^CMB ■ 



(4.6) 



(4.7) 



(4.8) 



Notice that Xcmb depends on the parameters {flmo, 
IfeOj h). The density parameter of radiation today is 

a-o = fi7o(l + 0.22717Vcff), (4.9) 



where Q-yo is the photon density parameter and N^s is the 
relativistic degrees of freedom. We adopt the standard 
values n^o = 2.469 x IQ-'^ h'^ and A^eff = 3.04 

The overall likelihood function'' is given by the product 
of the individual likelihoods according to: 



'Ctot(Pl,P2) = ^siPi) X /:/(pi,p2) 



(4.10) 



Zcmb 1048 (1 + Qmi24uj-" '^^^) (1 + giLj. 



0.0783a;^"°-2^V (l 

^b ji 



where gi 

0.560/ (l + 21.1a;i-«i), ujb = VLboK^ , and w,, 



39.5C.0-763), 52 



a-) , (4.4) 
VLmoh"^ {h 



^ Likelihoods are normalized to their maximum values. In the 
present analysis we always report la uncertainties on the fitted 
parameters. 
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ii„, n,„ n,„ 

FIG. 2: Left: Likelihood contours for equal to 2.30, 6.18 and 11.83, corresponding to la, 2a and 3a confidence levels, in the 
(fimo, 7) plane using a ACDM expansion model. Middle and Right: Here we show the corresponding contours in the case of f{R) 
models (/iCDM - middle panel and /2CDM -right panel). In all cases the red point corresponds to (r2mO,7) = (0.272,6/11). 
In this plot and in the ones that follow we have set the parameters that are not shown (eg 6) to their best fit values for the 
corresponding model (see Table HIT)) . 



where C / refers to the dynamical probe likehhood fit and 

^e{Pi) = CsNIa X £bAO X ^CMB ■ (4.11) 

The vectors pi, P2 contain the free parameters of the 
/(i?) model and depend on the model. In particular, 
the essential free parameters that enter in the theoretical 
expectation are: pi = {flmo^b) and P2 = (tojTi)- Note 
that in the case of the ACDM we have pi = r^mo- Also, 
in all cases we have set as — 0.8. 

Since likelihoods are defined as £ oc exp (— this 
translates into an addition for the joint Xtot function: 



with 



xL(Pi,P2) = xKpi) +X/(Pl,P2) 



xl(Pl) - XSNIa + XbAO + XCMB 



(4.12) 



(4.13) 



The minimization of the Xtot done 
Mathematica"'"'^ . Note that the X/ is given by 



Nf 

E 



/(^Sobsi^t) -/cts(P1,P2,Zj) 



(4.14) 

where cr^ is the observed growth rate uncertainty. 

To this end since N/ufu > 40 we will use, the relevant 
to our case, corrected Akaike information criterion [73| , 
defined, for the case of Gaussian errors, as: 



AIC = X 



2 

min 



2n 



fit 



(4.15) 



where N — Nexp + Nf, Ufn is the number of free 
parameters. A smaller value of AIC indicates a better 
model-data fit. However, small differences in AIC are 
not necessarily significant and therefore, in order to as- 
sess, the effectiveness of the different models in repro- 
ducing the data, one has to investigate the model pair 



difference AAIC= AICj^ - AIC^. The higher the value 
of |AAIC|, the higher the evidence against the model 
with higher value of AIC, with a difference |AAIC|^ 2 
indicating a positive such evidence and |AAIC|^ 6 in- 
dicating a strong such evidence, while a value <, 2 in- 
dicates consistency among the two comparison models. 
A numerical summary of the statistical analysis for the 
background expansion models as well as for the various 
7(2;) parametrizations is shown in Table Hill 



1. Constant growth index 

First of all we utilize the Tq parametrization (7 = 70, 
71 = 0: see section 3A). Therefore, the corresponding 
/iCDM and /2CDM statistical vectors P2 contain only 
three free parameters namely, P2 = (pi,7o,0) where 
Pi = {flmOib). Accordingly, if we consider the ACDM 
model then pi = f2mOj implying that the vector p2 in- 
cludes two free parameters. 

Our main results are listed in Table IIIIl where we 
quote the best fit parameters with the corresponding la 
uncertainties, for three different expansion models. In 
Figure [5] we present the la, 2a and 3a confidence lev- 
els in the (flmOTl) plane. It becomes evident that us- 
ing the most recent growth data-set together with the 
expansion cosmological data we can place strong con- 
straints on (5^mO,7)- lu all cases the best fit value 
^mo = 0.272 ± 0.003 is in a very good agreement with 
that provided by WMAP9-fSPf -f ACT (flrao = 0.272; 
Hinshaw et al. |74|). 

Concerning the ACDM expansion model (see the right 
panel of FiglJ) our growth index results are in a gree ment 
within la errors, to those of Samushia et al. [68| who 
found 7 = 0.584±0.112 and to those of [39| who obtained 
n,no = 0.273 ± 0.011 and 7 = 0.586 ± 0.080. However, 
our best-fit value 7 = 0.597 ± 0.046 is somewhat greater 
from the theoretically predicted value of 7a — 6/11 (see 
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0.6 



0.8 



0.7 
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0.3 
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ACDM 



0.4 



0.6 



0.8 



0.7 



0.6 



0.3 



0.2 
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ACDM 



0.0 
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0.8 



FIG. 3: Comparison of the observed and theoretical evolution of the growth rate fag{z) = F{z)as{z). Top: The blue, green 
and red lines correspond to the best-fit To, Fi and F2 parametrizations with the background expansion given by ACDM, while 
the black line to the exact solution of Eq. (|3.1|l for fas{z) for the ACDM model for n„o = 0.273 [z!]- Middle and Bottom: 
The blue, green and red lines correspond to the best-fit Fo, Fi and F2 parametrizations with the background expansion given 
by /iCDM and /2CDM respectively, while the black line to the exact solution of Eq. p.lf) for fas{z) for the ACDM model for 
n,„o = 0.272 Q. In all cases we utihze as = 0.8. 



lines in the right panel of Fig. [5]). Such a small dis- crepancy between the theoretical ACDM and observa- 
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TABLE III: Statistical results for the combined growth data (see Table I): The 1"* column indicates the expansion model, the 
2"^* column corresponds to 7(2) parametrizations appearing in section 3A. 3'^'* and 4'""' columns provide the QmO and b best 
values. The 5*'' and G**^ columns show the 70 and 71 best fit values. In all cases we used as = 0.8. The remaining columns 
present the goodness-of-fit statistics (xmim AIC and |AAIC| — |AICa — AlCf(ji)\). 



Exp. Model Param. Model Jlmo 




b 












/ -L 






AIC 


|AAIC| 


ACDM 


To 


0.272 


± 


0.003 








0.597 


± 


0.046 









574.227 


578.227 







El 


0.272 


± 


0.003 








0.567 


± 


0.066 


0.116 


± 


0.191 


573.861 


579.861 


1.634 




E2 


0.272 


± 


0.003 








0.561 


± 


0.068 


0.183 


± 


0.269 


573.767 


579.767 


1.540 


/iCDM-[21] 


To 


0.271 


± 


0.003 


0.111 


± 


0.140 


0.598 


± 


0.046 









573.855 


579.855 


1.628 


El 


0.271 


± 


0.003 


0.109 


± 


0.142 


0.573 


± 


0.068 


0.097 


± 


0.195 


573.633 


581.633 


3.406 




E2 


0.271 


± 


0.003 


0.109 


± 


0.142 


0.579 


± 


0.070 


0.101 


± 


0.275 


573.585 


581.585 


3.358 


/2CDM-[22] 


Eo 


0.272 


± 


0.005 


0.292 


± 


0.647 


0.594 


± 


0.047 









574.178 


580.178 


1.951 


El 


0.272 


± 


0.005 


0.150 


± 


1.355 


0.567 


± 


0.066 


0.113 


± 


0.199 


573.857 


581.857 


3.630 




E2 


0.272 


± 


0.005 


0.149 


± 


1.261 


0.561 


± 


0.068 


0.179 


± 


0.279 


573.765 


581.765 


3.538 



tionally fitted value of 7 has also been found by otlier au- 
thors. For example, Di Porto & Amendola [33| obtained 
7 = O.GOtoit Gong [131 measured 7 = 0.64^°;^^ while 
Nesseris & Perivolaropoulos J2g , found 7 = 0.67j;Q^7. 
Recently, Basilakos & Pouri [37| and Hudson & Turnbull 
[ssj l using a similar analysis found 7 — 0.602 ± 0.05 and 
7 = 0.619± 0.054 respectively. In this context, Samushia 
et al. obtained 7 = 0.65 ± 0.05. 

Concerning the /iCDM model (see the middle panel 
of Fig. ^ the best fit parameter is 70 = 0.598 ± 0.046 
with Xmin — 573.855. Also in the case of Starobin- 
sky's /2CDM model (see Fig. [2]) the best fit param- 
eter is 70 = 0.594 ± 0.047 with xL« - 574.178. In 
Figure [21 we plot the measured fcrsobsi^) ^i^h the esti- 
mated growth rate function, fag{z) = F{z)as{z). The 
value of AICa(~ 578.227) is smaller than the corre- 
sponding f{R) values which indicates that the ACDM 
model (7A — 0.597) appears now to fit better than the 
/iCDM and /2CDM gravity models the expansion and 
the growth data. The |AAIC| = |AICa - AIC/^_,(ij)| val- 
ues (ie., ~ 2) indicate that the growth data are consistent 
with the /iCDM and /2CDM gravity models with a con- 
stant growth index. 



2. Time varying growth index 

Now we concentrate on the 7(2) parametrizations, 
presented in section 3A. In this case the free param- 
eters of the models are p2 = (rimo,7o,7i) and p2 = 
(f}^o,^',7o,7i) for the ACDM and /iCDM, /2CDM ex- 
pansion models respectively. 

In FigureSKACDM model), Figure [5] (/iCDM model) 
and Figure [5] (/2CDM model) we present the results of 
our statistical analysis for the Fi (left panel) and r2 
(right panel) parametrizations in the (70,71) plane in 
which the corresponding contours are plotted for la, 2a 
and 3(7 confidence levels. The theoretical (70 771) values 
(see section 3 A) in the ACDM and f{R) expansion mod- 



els indicated by the colored dots. Overall, we find that 
the predicted ACDM, /iCDM and /2CDM (70,71) so- 
lutions of the ri_2 parametrizations are close to the la 
borders (Axi^r — 2.30; see green sectors in Figs. [H [S] and 
ED- 

Below we briefly discuss the main statistical results: 

(a) Fi parametrization: For the usual A cosmology the 
likehhood function peaks at 70 = 0.567 ± 0.066 and 71 = 
0.116±0.191 with xL« - 573.861. The latter results are 
in agreement with previous studies 26, 35-37]. In the 
case of the /iCDM and /2CDM gravity models we find 
that (70,71) = (0.573±0.068,0.097±0.195) with xL« - 
573.633 and (70,71) = (0.567 ± 0.066, 0.113 ± 0.199) and 
with Xmin — 573.857 respectively. 

(b) r2 parametrization: The best fit values are: (i) for 
ACDM we have 70 = 0.561 ± 0.068, 71 = 0.183 ± 0.269 
ixLin - 573.767), (u) in the case of /iCDM model we 
obtain 70 = 0.579 ± 0.070, 71 = 0.101 ± 0.275 (xL« - 
573.585) and (iii) for /2CDM gravity model we find 70 = 
0.561 ± 0.068, 71 = 0.179 ± 0.279 (xLn - 573.765). 

Finally, as we have already mentioned in Tabic IIIIl 
one may see a more compact presentation of our sta- 
tistical results. In Figure [7] we present the evolution of 
the growth index for various parametrizations. In the all 
three cases of the concordance A and /2CDM cosmolog- 
ical models (see upper and bottom panels of Fig[7]) the 
relative growth index difference of the various fitted 7(2;) 
models indicates that the ri_2 parametrizations have a 
very similar redshift dependence for z < 0.5, while the 
r2 parametrization shows very large such deviations for 
z > 0.5. Based on the /iCDM gravity model (middle 
panel of Fig[7]) we observe that the ri_2 parametriza- 
tions provide a similar evolution of the growth index. 

Finally, in Fig. [5] we show the evolution of Ge//(z) for 
the two f{R) models considered in the text, /iCDM (left) 
and /2CDM (right). The fines correspond to Fo (blue), 
Fi (green), and r2 (red). As it can be seen, both cases 
predict little evolution for Geff{z) at late times, around 
1.2% for /iCDM and 0.5% for /2CDM. Furthermore, 
while Geff{z) shows almost the same evolution for all 
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70 70 

FIG. 4: The ACDM expansion model: Likelihood contours (for Ax^ ~ — 21n£/£max equal to 2.30, 6.18 and 11.83, corresponding 
to la, 2a and 3a confidence levels) in the (70, 71). The left and right panels show the results based on the ri_2 parametrizations 
respectively. We also include the theoretical ACDM (70,71) values given in section 3, Ei = (6/11, 71(6/11, £7mO,i)/)) and 

E2 = (70,6/, 71(70,6/, f^mO, 6/))- 




FIG. 5: The Likelihood contours for the /iCDM expansion model (for more definitions see caption of Fig|4]). Here the 
colored dots correspond to the theoretical /iCDM (70,71) pair provided in section 3A, Ei — (6/11, 7i(6/ll, fimo,!)/)) and 

E2 = (70, b/, 71(70,6/, f^mO, 6/)). 



three parameterizations ro.1,2 of /iCDM, this is not the 
case for /2CDM where Fq differs significantly from the 
other two. 



5. CONCLUSIONS 

It is well known that the growth index 7 plays a key 
role in cosmological studies because it can be used as 
a useful tool in order to test Einstein's general relativ- 
ity on cosmological scales. In this article, we utilized 
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70 70 

FIG. 6: The Likelihood contours for the /2CDM expansion model (for more definitions see caption of Fig|4]). Here the 
colored dots correspond to the theoretical /2CDM (70,71) pair provided in section 3A, Ei — (6/11, 7i(6/ll, ilmo.f)/)) and 

E2 = (70,!)/, 71(70,6/, f^mO, 6/))- 



the recent growth rate data provided by the clustering, 
measured mainly from the PSCz, 2dF, VVDS, SDSS, 
6dF, 2MASS, BOSS and WiggleZ galaxy surveys, in or- 
der to constrain the growth index. In particular, per- 
forming simultaneously a likelihood analysis of the re- 
cent expansion data (SNIa, CMB shift parameter and 
BAOs) together with the growth rate of structure data, 
in order to determine the cosmological and the free pa- 
rameters of the 7(2) parametrizations and thus to sta- 
tistically quantify the ability of j{z) to represent the 
observations. We consider the following growth index 
parametrization 7(2;) = 70 + Jiuiz) [where y(z) = 0, 
y{z) = z and 1 — a{z)]. Overall, considering a ACDM ex- 
pansion model we found that the observed growth index 
is in agreement, within la errors, with the theoretically 
predicted value of 7a — 6/11. Additionally, based on the 
Akaike information criterion we shown that for any type 
of 7(2) the combined analysis of the cosmological (ex- 
pansion-|-growth) data can accommodate the Hu-Sawicky 
and Starobinsky f{R) gravity models for small values of 
the deviation parameter b. 

Numerical Analysis Files: The Mathcmat- 
ica and data files used for the numerical anal- 
ysis of this study may be downloaded from: 
http: / /leandros. physics. uoi.gr/fr-constraints/probes. htm 
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Appendix A: Useful formulae 

Here we provide the exact expressions of the coeffi- 
cients 6Hf{N) for the HS and Starobinsky models for 
n = l. In all cases H\{N) is given by Eq. 07T2\ . For the 
HS model the first two terms are: 




FIG. 7: Top: The evolution of the growth index 7(2) — ^ for the ACDM model. The lines correspond to To (blue), Fi 
(green), and r2 (red). Middle and Bottom: The evolution of the growth index for the /iCDM and /2CDM cosmological models 
respectively. The lines correspond to Fo (blue), Fi (green), and F2 (red). The dashed lines correspond to the la errors. 
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FIG. 8: The evolution of the Geff{z) for the two f{R) models considered in the text, /iCDM (top) and /2CDM (bottom), for 
all three growth rate parametrizations To (left), Fi (middle), and F2 (right). 



6H^{N) ^ (-1 + n^o + ^ro) (6gA(^)' + 4.H'^{Nf + HAjN) {IbH'^jN) + 2H'j[{N))) 

Hi 2HA{N){2HA{N) + H'j^{N)f 



(A.1) 



{H^{-i + n„ + ^lr)\l2mK{Nf - s2hI{-i + + nr)H'j,{Nf + 32irA(iV)^(25ir;(7V) 



+3H'l{N)) - 2H^{-1 + + nr)HAiN)H'j,{N)^il39H'j^iN) + 22H'i{N)) 
+16HA{Nf{9H^{-l + + aO + 89H[^{Nf + 12H'a{N)H'j;^{N)) 

+HA{N)^H'jy{N)^{-7i9H^{-l + a„ + nr)H'jy{Nf + 9H'a{N)^ - A8H^i-l + a„ + nr)H'ji{Nf 
-4H^{-1 + n,n + nr)H'AiN)i74H'^iN) - SHa^^^N))) + SHAiNf (lUH^i-l + + ^MiN) 
+U6H'A{Nf + 18H'AiNfH'j[{N) + H^i-l + 0„ + f2,)(15if;^(iV) - 67?a(3)(7V) - Ha^^\N))) 
+4Ha{N)^{5A0HI{-1 + flm + flr)H'AiNf + 12AH'A{Nf + 12H'A{Nf Hl{N) 
+3Hl{-l + 0„ + nr)H'l{N){nH'i{N) + AHa^^\N)) + 2Hl{-l + + nr)H'A{N){l29H'l{N) 
+12iJA(3)(iV) - Ha^'^\N))) - 2HA{Nf{MHl{-l + a„ + nr)H'A{Nf - 53H'a{N)'' 
-3H'a{N)^H'1{N) + 21ij2(_l + 0„ + nr)H';{Nf + 3H^{-1 + + nr)H'A{N)H'^{N) x 
{ilH^liN) - iHA^'^HN)) + Hl{-1 + r!„ + ^lr)H'A{Nf {211 H'i{N) - A2Ha^^\N) + Ha^^\N))))) 
l{mA{N)\2HA{N) + H'A{N)f) (A.2) 



while for the Starobinsky model the first two terms are: 



5m{N) ^ (-1 + ^mo + »ro) ' {8HA{Nf + 9H'j,{Nf + Ha{N) (34gj,(iV) + QH'i{N))) 

Hi 4HA{N)^2HA{N)+H'^{N))^ ^ ""^^ 
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^2 = (-1 + ^mo + flro) ^{192HA{Nf - 486ff2 (_i + + fi^o) H'j^{N)'' + 320HAiNyi6H'j,{N) - 

+H'ji{N)) ~ 18H^ (-1 + n^o + ^ro) Ha{N)H'^{N)^ {263H'jy{N) + ^8H';{N)) + 16HA{Nf{32 (-1 + il^o + ^ro) + 

5H'j,{N) {47H'j^iN) + SH'j^iN))) + HA{NfK{Nf{-15562H^ (-1 + QmO + i^ro) H'^{Nf + 25H'j,{Nf - 

972ij2 (_i + + fi^p) H';{N f - l2Hl (-1 + il^o + ^^ro) K{N) (b2,2H';^{N) - lbHA^''\N)^ ) + 

AHa{N)\AA22HI (-1 + f]„o + ^ro) KiNf + 3A5H'f,{Nf + 40H'^{Nf H';^{N) + 

12El (-1 + 17^0 + f^ro) {9H'^{N) + 2ifA'''(A^)) + ISH^ (-1 + f^mo + ^^ro) K{N) x 

{l75H';{N) + 20Ha^^\N) - iJA^^'w)) + 2HA{Nf{-7000H^ (-1 + 0„o + ^^ro) KiNf + 

U8H'^{Nf + 10H'f,iN)^H';{N) - S24H^ (-1 + n^o + ^ro) K^Nf - 36H^ (-1 + + f^ro) H'^{N)H';{N) x 

(6iH';{N) - 4FA(3)(iV)) - 9H^ (-1 + n^o + ^ro) H'^{Nf {&mH'i{N) - &&Ha^^Hn) + Ha^^Hn))) + 

8HA{Nf{82m^ (-1 + + ^ra) + ^OOH'j,{Nf + 60H'j,{N)^H';{N) + 3H^ (-1 + Qrno + ^ro) x 

(29H'^{N) - 3 {6Ha^^\N) + Ha^^\N)))))/ {l6HA{Nf {2Ha{N) + H'^{N)) 1°) (A.4) 
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